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^ ■ We use ISAJET to perform a detailed study of the multilepton signals 

. expected from cascade decays of supersymmetric particle produced at the 

I CERN LHC. Our analysis is performed within the framework of the minimal 

QhI supergravity model with gauge coupling unification and radiative electroweak 

Q^' symmetry breaking. We delineate the regions of parameter space where jets 

^ ■ plus missing energy plus 1, 2 (opposite sign and same-sign dileptons), and 3 

• • . isolated lepton events should be visible above standard model backgrounds. 

I We find that if any I^j, signal at the LHC is to be attributed to gluino and/or 

^ ' squark production, and if rUg ^ 1 TeV, then several of these signals must be 

5— ( ' _i 

' simultaneously observable. Furthermore, assuming 10 fb of integrated lumi- 

nosity, we find that the reach for supersymmetry in the l£+jets + ^j. channel 
extends to nig ~ 2300 (1600) GeV for rUg ~ nig {nig ~ l.bnig), and exceeds 
the corresponding reach in the Oi + I^j, channel. We show that measurements 
of the various topological cross sections, jet and 5-hadron multiplicities in 
these events, together with the charge asymmetry for single lepton and same- 
sign dilepton events, and flavor asymmetry for opposite sign dilepton events, 
serve to narrow the allowed range of underlying SUGRA parameter values. 
We also delineate parameter regions where signals with clean isolated dilep- 
ton (from slepton production) and trilepton events (from chargino/neutralino 
production) are visible at the LHC, and examine the extent to which these 
signals can be separated from other SUSY sources. 
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I. INTRODUCTION AND MOTIVATION 



The search for supersymmetric particles is now an integral part of all current, as well 
as future, experimental programs at high energy colliders. Aside from the many attractive 
features of supersymmetry (SUSY), the impetus for these searches comes from the fact 
that weak scale SUSY which is introduced to ameliorate the fine tuning problem of 
the Standard Model (SM), requires that the supersymmetric partners of SM particles must 
be accessible to experiments that probe the TeV energy scale. Thus, while experiments at 
LEP2 and at the Tevatron (or its upgrades) may well discover sparticles, a definitive search 
for supersymmetry can only be performed at supercolliders such as the CERN Large 
Hadron Collider (LHC) or, at electron-positron linear colliders with ^/s = 500 — 1500 GeV. 
There is general agreement, based on detailed studies of SUSY signals both within the more 
general Minimal Supersymmetric Model (MSSM) framework as well as within the very 
attractive and economic supergravity (SUGRA) GUT framework ||^-|TT|, that weak scale 
SUSY will not evade detection at these facihties 

The natural question then is: if we do see signals for new physics, can we unravel 
their origin, and trace them to the production of supersymmetric particles? At electron- 
positron colliders where the cleanliness of the interaction environment allows for [|I4|,|8|,|9|JT5 



the precision measurement of at least some of the properties of these particles (mass, spin, 
decay patterns,...), this may be a straightforward exercise, especially if the machine energy 
is increased incrementally, so that it is possible to focus on just one new signal at a time. 
At the LHC, however, the situation is much more complicated, not only because of the 
messier environment, but also because all new particles which are kinematically accessible 
will simultaneously contribute to the signal: we will thus have the additional task of sorting 
the supersymmetric signals from one another in order to discover the nature of the new 
physics. 

Some progress has already been made on the issue of identifying the sparticle production 
processes that give rise to SUSY signals at the LHC. For instance, it has been shown 
that, with suitable cuts, the clean 3i + J^j, signal from the production of charginos and 
neutralinos via the reaction pp W1Z2 ivZ\ + I'l' Z\ can not only be separated from 
SM backgrounds, but also, that it can be isolated from other SUSY sources. An observation 
of a signal in this channel would, therefore, unambiguously point to W\'Z2 production as its 
source, at least within the SUSY framework. It is, however, not always possible to devise 
cuts to isolate a single source of SUSY events. A detailed study of the signal characteristics 
may then help to identify the sparticles producing the signal. In a previous study [1TT| , 



hereafter referred to as Paper I, we examined the reach of the LHC in the multijet plus 
l^j, channel and studied what information could be obtained by a detailed study of this 
sample. Assuming as usual that squarks cannot be much lighter than gluinos, we showed 
that if gluinos are lighter than about 750 GeV, their mass could be extracted to 15 - 25% by 
reconstructing multijet masses in opposite detector hemispheres. Furthermore, by measuring 
the mean jet multiplicity, (n^), which is observably larger if squarks are much heavier than 
gluinos, it should be possible to distinguish the ~ m-g case from the one where squarks 
are substantially heavier than gluinos. While the gluino mass would determine mi/2, 
will at least enable us to decide whether mo is small or comparable to mi/2, or much larger. 
Are there other ways by which we can tell what is being produced at the LHC? Also, 
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is it possible to test whether the minimal SUGRA framework adopted in Paper I (as well 
as in many recent phenomenological analyses @J^,^^23| of SUSY) can consistently account 



for all the observed signals, or whether some of the underlying assumptions about the sym- 
metries of physics at the ultra-high scale need to be modified? It has already been shown 
that the precision measurements that are possible in the clean environment of e^e~ 
collisions will allow experimentalists to perform incisive tests of the SUGRA framework at 
future Linear Colliders. While it is not possible to perform similar measurements at hadron 
colliders, the big advantage of the LHC over 500 GeV Linear Colliders is that many more 
sparticle production processes should be kinematically accessible, resulting in a large number 
of potential observables. Since the minimal SUGRA GUT model with radiative breaking 
of electroweak symmetry is completely fixed by just four SUSY parameters: mo and mi/2, 
the universal scalar mass and gaugino masses at the high scale Mx ~ Mqut, the SUSY 
breaking universal trilinear coupling Aq, and the parameter tan/3 along with sgn /i, the con- 
sistency of the framework can be tested by verifying that the rates and distributions in all 
the observed channels can be accommodated by a single choice of model parameters. Even 
more ambitiously, one could ask whether it would be possible to determine the underlying 
parameters from the observed signals, and we report the results of our preliminary attempt 
to do so in this paper. 

We stress here that we do not mean to imply that the SUGRA framework is the uniquely 
correct one. Indeed the sensitivity to the details of its predictions should be examined, 
particularly when studying the reach of future facilities. Nonetheless, it is an economic, 
attractive and predictive framework, and it can be used as a guide for sparticle masses and 
mixing patterns. Such a framework is needed since without assuming anything other than 
the weak scale symmetries, there are far too many parameters, making phenomenological 
analyses intractable. 

It should be clear from the preceding discussion that a study of all possible signals as 
a function of SUGRA parameters is a first step toward testing the model framework at the 
LHC. Of course, it is equally important to quantify the reach of the LHC in each of these 
channels which include, 

• the non-leptonic I^j, + jets channel studied in Paper I, 

• the li + jets + I^j, channel, 

• the opposite-sign (OS) dilepton + jets + I^j, channel, 

• the same-sign (SS) dilepton + jets + channel, 

• the multilepton + jets + I^j, channel, with ne > 3. 

These signal channels ought to originate mainly from squark and gluino pair production, 
followed by their cascade decays. In addition, there are also clean [i.e. free from central jet 
activity) channels with 

• dilepton plus events, and 

• trilepton and n£ > 4 lepton events. 
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mainly from the pair production of sleptons [H,^] as well as from WiWi, W1Z2 |l26|JT 
and Z2Z2 ^\ production processes. Of course, after cuts, several sparticle production 
mechanisms can contribute to each channel, so that it is necessary to simultaneously generate 
the production of all sparticles in order to obtain an accurate assessment of the expected 
signals. 

In this paper we continue the study of SUSY signals within the SUGRA framework that 
we began in Paper I. We use ISAJET 7.14 p8[ to compute p9[ the signal cross sections after 
cuts designed to separate the SUSY signals from SM backgrounds, and wherever possible, 
also to separate SUSY sources from one another, in each of the leptonic channels listed 
above. We show these in the mo - mi/2 plane, which provides a convenient way to display 
the signals from different sparticle production processes. For other parameters, our canonical 
choices are Aq = 0, tan/3 = 2 and 10 and we adopt both signs of /i. To orient the reader 
with various sparticles masses derived from the SUGRA framework, we show contours of 
squark and gluino masses in Fig. |l|, and of slepton and chargino masses in Fig. ^ for 
(a) tan/3 = 2,/i < 0, (b) tan (3 = 2,/i > 0, (c) tan/? = 10,/i < 0, and (d) tan/3 = 10,/i > 0. 
We remind the reader of the approximate relationship m^^ ~ "^#1 — '^^Zi ^^^^ usually holds 
because |//| tends to be large within this framework. In Fig. |^, as well as in many subsequent 
figures, the regions shaded by bricks (hatches) are excluded by theoretical (experimental) 
constraints as discussed in Paper I. The gluino mass contours in Fig. |l| are not exactly 
horizontal because of the difference I^O] between the running and physical {i.e. pole) gluino 
mass. In Fig. H, we also show the region where the "spoiler" decay modes Z2 —>■ ZiHi or 
Z2 ZiZ are kinematically accessible (above the dotted contours); in this region, leptonic 
decays of the Z2 are either very suppressed, or have additional backgrounds from SM Z 
boson production. 

We map out the regions of parameter space where these signals are observable at the 
LHC, and compare this with the region that can be probed via the Ifrp channel [11,3| as 
delineated in Paper I. On the issue of the LHC reach, our main new result is that the l^ 
channel provides the greatest reach for supersymmetry. However, the observation of signals 
in several channels is important, since it can help to identify SUSY as the unique source 
of new physics. We study jet and 5-hadron multiplicity distributions, as well as charge 
asymmetry distributions in the single lepton and SS dilepton + jets + I^j, channels, and 
dilepton flavor asymmetry in the OS dilepton + jets + I^j. channel as these can provide 
information about the cascade decay chains of gluinos and squarks [0]. We also identify 
regions of parameter space where the clean trilepton and the clean OS, same-flavor dilepton 
signals are observable. While these regions form a subset of the region where SUSY may be 
probed via the multijet channels, an observation of these signals will be important because 
they will signal W1Z2 and slepton production, respectively; i.e. with suitable cuts described 
below, there is limited contamination from other SUSY sources. 

The rest of this paper is organized as follows. In the next Section, we briefly discuss some 
computational details. Sec. 3 and Sec. 4 focus on the multilepton plus multijet and clean 
multilepton channels, respectively. We present a comparative analysis of the reach in various 
channels in Sec. 5, and also consolidate the information about the underlying SUGRA model 
parameters that might be obtained by studying distributions in these various channels. 
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II. EVENT SIMULATION 



We work within the framework of the minimal SUGRA model and use IS A JET 7.14 
to simulate the various leptonic signals for SUSY listed above. The implementation of the 
SUGRA framework into ISAJET has been described elsewhere [pO| , pT| and will not be re- 
peated here. We generate a// lowest order 2 — > 2 SUSY subprocesses in our simulation of the 
multilepton plus multijet signals (except for s-channel Higgs boson mediated subprocesses). 
However, for the simulation of the clean multilepton signals, we have generated only slepton 
and chargino/neutralino events, since gluino and squark decays will very seldom yield final 
states without central jet activity . 

For detector simulation at the LHC, we use the toy calorimeter simulation package ISA- 
PLT. We simulate calorimetry covering —5 < rj < 5 with cell size At] x A0 = 0.05 x 0.05. 
We take the hadronic energy resolution to be 50%/\/E © 0.03 for Ir^l < 3, where © denotes 
addition in quadrature, and to be 100%/\/E © 0.07 for 3 < |?7| < 5, to model the effective 
Pt resolution of the forward calorimeter including the effects of shower spreading, which is 
otherwise neglected. We take electromagnetic resolution to be 10%/ © 0.01. Although 



we have included these resolutions, which are typical of ATLAS and CMS [0, we have 
made no attempt to estimate the effects of cracks, edges, and other problem regions. Much 
more detailed detector simulations are needed to understand the effects of such regions and 
of the resulting non-Gaussian tails, particularly on the resolution. 

Jets are found using fixed cones of size R = ^/Arf~+~A(p = 0.7 using the ISAJET routine 
GET JET. Clusters with Ex > 100 GeV and |?7(jet)| < 3 are labeled as jets. However, for 
the purpose of jet- veto only, clusters with Et > 25 GeV and |?7(jet)| < 3 are regarded as 
jets. Muons and electrons are classified as isolated if they have pt > 10 GeV, \ri{£)\ < 2.5, 
and the visible activity within a cone of i? = 0.3 about the lepton direction is less than 
ET(cone) = 5 GeV. 

We assume an integrated luminosity of 10fb~^, corresponding to 10'^^cm~^s~^ for one 
year. Hence we feel justified in neglecting the effects of pileup. We presume it would 
be possible to use the maximum LHC luminosity, lO^'' cm~^s~^, to search for gluinos and 
squarks with masses ^1-2 TeV. 



III. MULTILEPTON PLUS MULTIJET SIGNALS FOR SUPERSYMMETRY 
A. Classification of Signals and Event Selection 

For rUg, nig < 1 TeV, gg, gq and qq production is the dominant source of SUSY events 
at the LHC. These production mechanisms, together with g and q cascade decays, naturally 
lead to events with nleptons + mjets + I^j,, where typically n = - 4 and m > 2. These 
event topologies may also arise from the production of gluinos and squarks in association 
with a chargino or a neutralino. In addition, direct production of charginos, neutralinos and 
sleptons followed by cascade decays to Wi or Zj can lead to similar events. 

Although in our simulation we generate all SUSY processes using ISAJET, our cuts are 
designed to selectively pick out gluino and squark events, whose characteristics are high 
transverse momentum jets and large missing transverse energy. Furthermore, the pt of the 
primary jets from gluinos, as well as the are expected to scale with rrig. In contrast. 
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the momenta of leptons, produced far down in the cascade decay chain from chargino and 
neutrahno daughters, will not scale in energy the same way as jets and Ifj, which can be 
produced in the first step of the cascade decay. Thus, following Paper I, for the multilepton 
plus multijet signals for SUSY, we vary the missing-energy and jet Et cuts using a parameter 
but fix the lepton cuts: 

• jet multiplicity, n^et > 2 (with -Erjet > 100 GeV), 

• transverse sphericity St > 0.2, 

• EtUi), Et{]2) > and > E^. 

We classify the events by the multiplicity of isolated leptons, and in the case of dilep- 
ton events, we also distinguish between the OS and the SS sample as these could have 
substantially different origins. For the leptons we require, 

• prii) > 20 GeV (£ = e or /i) and Mt(£, > 100 GeV for the li signal, and 

• Pt(^i, ^2) > 20 GeV for n = 2, 3, . . . lepton signals. We do not impose any pt(^) > E!j. 
requirement on the leptons for reasons explained above. 



B. Calculation of Backgrounds 

SM processes, particularly those involving the production of heavy particles like the W 
and Z bosons, or the top quarks, can mimic the leptonic signals listed above. We have used 
ISAJET to evaluate the following SM backgrounds to these signals: 

1. ti production, where the leptonic decays of the tops can give up to two isolated leptons; 
for n > 2 the additional lepton may come from a 6 or c decay or from the fragmentation 
of additional jets in the event, where the lepton is accidentally isolated; 

2. W and Z boson + jet production, where additional jets and/or leptons come from 
parton showering; 

3. WW, WZ and ZZ production, where additional jets can again arise from QCD radi- 
ation; 

4. QCD jet production, where leptons can arise from decays of heavy flavors produced 
directly or via gluon splitting. 

ISAJET includes higher order QCD and electroweak effects in the branching approxima- 
tion: i.e. it includes quark and gluon as well as weak vector boson radiation, using exact 
kinematics but only coUinear dynamics. Thus, extra leptons can arise in any of the above 
hard scattering subprocesses additionally from, for instance, gluon splitting to top or bottom 
quark pairs, followed by their subsequent decays, or by W and Z boson radiation. 

The Efj. dependence of these background cross sections, obtained using CTEQ2L par- 
ton distributions |Q, is displayed in Fig. || for (a) 1£ + jets events, (b) OS dilepton -|- jets 
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events, (c) SS dilepton + jets events, and (b) 3£ + jets events. We see that W + jets pro- 
duction is generally the largest background, except in the OS dilepton channel where the 
tt background dominates for modest values of E'f^. Gauge boson pair production and QCD 
background sources are essentially negligible, compared to backgrounds from W, Z and 
ti production. The multilepton background from gauge boson pair production is strongly 
suppressed, presumably because of the requirement of two additional hard jets as well as 

The wiggles in the curves in Fig. ^ are a reflection of the statistical fluctuations in our 
simulation. We see that for modest values of the fluctuation in the biggest backgrounds 
are under control. We will see later that in order to extract the reach we use E^j. = 200 GeV 
in all but the li channel for which we use E^ = 400 GeV. We have checked that for these 
ranges we typically obtain at least several tens of events passing the cuts in our simulation, 
so that the statistical errors on the relevant background estimates are in control. It should, 
of course, be remembered that our background calculations are probably correct only to a 
factor ~ 2-3 due to the inherent uncertainties associated with leading-log QCD, the parton 
shower approximation, our idealistic detector simulation, etc. 

In order to enable the reader to assess this calculation, we have shown the details of the 
background calculation in the various multilepton channels in Table I for one value of E^. 
Since only a tiny fraction of the events generated pass the cuts, it is necessary to generate 
events in several ranges of hard scattering px {p^^) for each SM process, and then combine 



these to obtain the background cross section from each of these sources |34]. The results of 



our computation for E^ = 200 GeV are shown in Table I, for the li, OS, SS and 3i signals. In 
those p^''^ bins where we obtain no events, the bound shown corresponds to the cross section 
corresponding to the one-event level. We see that for the tt and W or Z backgrounds 
which are the largest contributors to the background cross section, the main contribution 
indeed comes from the intermediate values of pff^, ensuring that we do have a reasonable 
estimate for the cross section. We have also checked that for the major contributors to the 
background, we have ten to several hundred events passing the cuts in our simulation, so 
that our estimates should be reliable to a few tens of percent, and frequently much better, 
as far as statistical errors are concerned. Finally, we see that the QCD background to the 
li cross section is clearly small; while we typically obtain only a bound on this from our 
simulation, it is reasonable to expect that this will not be a substantial background in the 
multilepton channels. 

We have also attempted to estimate the Ai background with ISAJET. Such events, how- 
ever, form an extremely tiny fraction of the total cross section so that a reliable simulation 
of these would require lengthy computer runs. Our simulation in which just a handful 
of events pass the cuts in each of the iy-|-jets, Z-l-jets and the ti channels yields a cross 
section a{Ai) = 0.04 fb for this background for E^ = 100 GeV, which falls to 0.002 fb 
for E^ = 200 GeV. Even allowing for uncertainties in our estimates, we see that the SM 
background is essentially negligible. For reasons of brevity and because these signals are 
observable only for limited ranges of parameters, we have not shown these cross sections in 
the figures or in Table I. 
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C. SUSY Multilepton plus Multijet Signals at the LHC 



Our next goal is to evaluate the various SUSY ra-lepton + m-jets + Ifj, signals expected 
from supersymmetry at the LHC, and compare against background expectations. Toward 
this end, we show the signal cross sections along with the total SM background as a function 
of in Fig. ^ for (a) 1£ + jets events, {h) OS dilepton + jets events, (c) SS dilepton + jets 
events, and (6) 3£+jets events. Our total signal and background cross sections are evaluated, 
as usual, at leading-log level, and so are uncertain to about a factor of 2; next-to-leading log 



gluino and squark cross sections can be found in Ref. |25|- We have illustrated the signal 



for the same six choices of SUSY parameters as in Ref. UTTI]; we take Aq = 0, tan/? = 2, 



mt = 


170 GeV and 




1. 


mo 


= mi/2 = 


100 GeV, for which rrig = 


2. 


mo 


= 4mi/2 = 


= 400 GeV, for which rrig = 


3. 


mo 


= rni/2 = 


300 GeV, for which rrig = 


4. 


mo 


= 4mi/2 = 


= 1200 GeV, for which rrig 


5. 


mo 


= "^1/2 = 


600 GeV, for which rrig = 


6. 


mo 


= 4mi/2 = 


= 2000 GeV, for which rrig 



1300 GeV and m,- = 2200 GeV. 

From Fig. ^ it is relatively obvious how should be chosen to search for SUSY in the 
multilepton plus multijet channels: if gluinos are relatively light (cases 1 and 2), E^ ~ 100 
-150 GeV suffices to obtain a large signal to background ratio and a large event rate in all 
the channels. For the cases with heavier gluinos and squarks (cases 3-6), a larger value of 
E^ is necessary, though it should not be chosen too large as to cut out all the signal. For 
instance, Efj. ~ 200 GeV should yield an observable signal, with a signal to background ratio 
larger than unity in all but the OS dilepton channel. The maximal reach may be anticipated 
to occur in the li channel — for cases 5 and 6, with E^ = 400 GeV, we expect ~ 20 - 
100 events (versus a background of just about three or four events) after a year of LHC 
operation at its "low" luminosity of 10fb~^/yr. 

We next examine in detail each of the multilepton plus multijet topolog function 
of SUGRA parameters. 



1. Single Lepton Events 

We begin by showing, in the mo - mi/2 plane, cross section contours for the li signal 
after the cuts discussed above for Aq = and (a) tan/? = 2, < 0, (6) tan/? = 2, > 0, 
(c) tan/? = 10, /i < 0, and (d) tan/? =10, /i > in Fig. ^. We have shown the results 
for Ej, = 100 GeV (solid) for which the total SM background from Fig. ^ is ~ 1300 fb, and 
also for E^ = 400 GeV (dotted), for which the background is very tiny at about 0.5 fb. 
For an integrated luminosity of 10fb~^ the corresponding 5a limits are, 57 fb and 1.1 fb, 
respectively. 
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To obtain these contours (as well as the corresponding contours for the OS, SS and 3£ 
signals discussed below) we have first computed the signal cross section for each point on a 
100 GeV X 100 GeV lattice in the mg - '"^1/2 plane, for points which do not fall inside the 
excluded shaded regions. The contours are then obtained via interpolation. We have cut off 
the contours near the boundaries of the shaded regions where the sampling is poorer and 
the interpolation not as reliable. 

The curves shown are for cross sections of 1, 2, 4, 8, . . .fb (only every other solid curve 
is labeled) in each of these cases. We have also checked that, even for the very hard 
cut, there are sufficiently many events in our simulation to yield reliable estimates of the 
cross sections: for the Efj. = 400 GeV case, the efficiency for SUSY events to pass the cuts 
becomes very small unless sparticles are rather heavy, so that for moderate mi/2 values very 
lengthy computer runs would be necessary to compute the cross section. For this reason, 
and because very hard cuts are necessary only for the largest gluino and squark masses, we 
have shown only the first three dotted curves, corresponding to cross sections after cuts of 
1, 2 and 4 fb. 

We see from Fig. ^ that with = 100 GeV, the 5a reach (the 64 fb contour is closest 
to the 57 fb 5a limit) in the li channel extends to mi/2 ~ 600 GeV for small values of mo 
(corresponding to nig ~ nig ~ 1.5 TeV), or to mi/2 ~ 400 GeV {rrig ~ 1 TeV) if squarks 
are heavy. Regions below this 5a contour all have larger signal cross sections, so that we 
found no "holes" of non-observability below the 5a limit. Notice, however, that the signal 
to background ratio at the 5a limit is just less than 5%: if we require this ratio to exceed 
25%, the corresponding reach is between the 256 and 512 fb contours. To probe values of 
mi/2 ~ 250 — 300 GeV, it is best to choose larger values of Efj. to obtain better statistical 
significance as well as higher signal to background ratio. The maximal reach in the li channel 
can be obtained by using a hard Efj. cut which eliminates essentially all the background but 
still retains the signal at an observable level. The highest of the dotted contours (1 fb) 
is very close to the 5a limit for E^ = 400 GeV; in this case, signal/background ~ 1, and 
^1/2 ~ 700 GeV (1000 GeV) can be probed in the large (small) mo region. This corresponds 
to a reach in ~ 1700 GeV (2300 GeV). Thus, we note that the reach in this channel 
appears to substantially exceed the corresponding reach in the canonical multijet + I^j, 

(no isolated lepton) channel. 

How well can one determine the SUSY parameters by studying the 1£ + jet + I^j, signal? 
Measurement of the total rate for such events would localize, within errors, a position along 
one of the total cross section contours of Fig. ^. These contours vary strongly with mi/2, 
but less strongly with mg. The cross sections are roughly the same in all the four frames. 
While this means that the 1£ signal rate yields no information about tan /? or sgn /x, it also 
means that a measurement — and calculation — of the cross section to within a factor of 2 
would indeed tell us on which contour we are within about ±50 GeV. While this does not 
accurately pin mi/2 because the contours are not quite horizontal, one would still be able 
to obtain a reasonable estimate of mi/2- The range of mi/2 thus obtained can be checked 
for consistency with the gluino mass that might be extracted [|r^ from the I^j. channel: in 
fact, a similar measurement ought to be possible in the 1£ channel. 

Determination of mo is more difficult and will probably require a simultaneous study of 
several signals and their distributions. We note, however, that for very small values of mo, 
because of the enhancement of the leptonic decays of Wi, and frequently also of Z2, the 
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lepton plus multijet cross sections are large. In contrast, the Oi plus multijet cross sections 
(for a fixed value of mi/2) actually reduce 0] as rriQ because of the lepton veto: thus, a 
measurement of the ratio of non-leptonic to multileptonic multijet cross sections could yield 
information on whether we are in the small mo region, particularly in the region where two 
body decays of Wi and/or Z2 into real sleptons are kinematically allowed (this region has 
been delineated in Fig. |15| below). 

Next, we turn to an examination of the single lepton charge asymmetry, which could 
provide an additional handle on position in parameter space. Since LHC is a pp collider, 
there is a preponderance of valence w-quarks in the initial scattering state, which can lead to 
a large proportion of u squarks being produced in the final state if squarks are moderately 
heavy and gluinos not too light (otherwise, sea parton annihilation may be the dominant 
source of sparticles). The -u-squarks frequently decay to W^, and one is led to expect more 
i~^s being produced than i~s. For large ttiq compared to mi/2, gg production should be 
dominant, which leads to equal production of and £~ in cascade decays. Likewise, as 
we already noted, if mi/2 is small, then sparticle production at the LHC is dominated by 
gluon fusion and sea-quark annihilation, which also leads to equal and i~ production. In 
contrast, for larger values of mi/2 and not too large mo, the squark production via valence 
quarks can dominate, and lead to the lepton charge asymmetry. To illustrate this, we show 
in Fig. ^ the single lepton charge asymmetry 

_ N{i+) - N{i-) 
' ~ N{i+) + N{£-) ' 

versus mo, for (a) mi/2 = 100 GeV, using = 100 GeV, (b) mi/2 = 200 GeV, using 
= 200 GeV, (c) my 2 = 400 GeV, using E^ = 400 GeV and (d) mi/2 = 500 GeV, using 
E^ = 400 GeV. In all frames, Aq = and tan/3 = 2. For frames (a - d), we take /i < 0; 
frames {e)-{h) are the same except that /i > 0. The horizontal dashed line is at Ac = 0. The 
rather small SM background (see Fig. |^) has not been included in these figures. We indeed 
see that for the small mi/2 cases of frames (a) and (e), the asymmetry is consistent with 
zero (by choosing a larger E^ value, it may be possible to enhance the valence contribution 
and so obtain an asymmetry even in this case.). As we move up in mi/2 values, a significant 
positive charge asymmetry develops, especially for small values of mo, refiecting the relative 
contribution of M-squarks versus d-squarks, or other squark fiavors or gluinos. Knowledge 
of mi/2 may thus be combined with the measurement of Ac to roughly localize mo — we 
would at least learn whether mo ^ mi/2 or whether mo ^ mi/2. 

In order to explore other strategies for the determination of mo, in Fig. ^we have shown 
the mean jet multiplicity {{uj)) as a function of mo for the same eight cases (a) - (h) as 
in Fig. ^. We see that for a fixed value of mi/2 (which can be determined from other 
considerations), {rij) clearly increases with mo. The underlying physics is exactly the same 
||Tl| as for the sample: for small mo, squark production is a significant source of J^rp 
events, and because qr frequently directly decay via qn — > qZi, the mean jet multiplicity 
is reduced. The mean jet multiplicity is essentially independent of the sign of /x. We see, 
however, that it can increase by as much as a whole unit as mo varies between 100 GeV 
and 1 TeV. We also see that the precision with which mo can be determined depends on the 
values of other SUSY parameters. Finally, we note that although {rij) changes only by about 
30% as mo is varied over the whole range in the figure, this could mean a significant increase 
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in the cross section for high multiphcity (say nj > 4 or 5) relative to nj = 2, so that ratios 
of cross sections with different jet multiphcities could yield a more sensitive measure of rriQ. 
Devising the optimal measure for localizing ttlq would require a detailed study beyond the 
scope of the present analysis. 

The multiplicity of tagged S-hadrons may also yield information about the underlying 
parameters. Towards this end, in Fig. ^ we have plotted the mean multiplicity {{tib)) of 
tagged i?-hadrons in the 1£ SUSY sample for the same cases as in Fig. assuming that 
a i?-hadron with px > 20 GeV and |?7b| < 2 is tagged with an efficiency of 40%. We see 
from the figure that (ub) varies between 0.2 to 1.3 over the parameter range shown. For 
the light gluinos cases with mi/2 = 100 GeV in frames (a) and (e), (n^) is small and shows 
little variation with mo, except around mo = 200 GeV where the decays g bb dominate 
other squark decays. For large mo, the gluino decays via the three body modes, except 
that the decays to tops (which can potentially be enhanced) are kinematically suppressed. 
For heavier gluinos {i.e. larger values of mi/2) there are two important differences. First, 
the spoiler decay Z2 —>■ ZiH^, which is a source of 5's, may be kinematically accessible: 
this leads to an increase in {ub) which is roughly independent of mo, except for the very 
small mo region where neutralino decays to sleptons are also accessible. Second, gluinos 
are heavy enough to decay to t quarks. Thus, when mo is very small, gluinos dominantly 
decay via g ^ qq into all flavors. As mo is increased, g tit or bib may be kinematically 
allowed, while the g ^ qq modes are closed. For even larger values of mo, three body 
decays to third generation quarks can be enhanced due to propagator and large Yukawa 
coupling effects 0,^, leading to an increase in (ub). Finally, we remark that (n^) does 



not serve to discriminate between the two signs of /i. We caution the reader that {tib) may 
potentially be sensitive to variations in Aq, since these may alter the masses and mixings of 
third generation sfermions. Thus some care must be exercised when attempting to extract 
mo from a measurement of the i?-hadron multiplicity. 



2. Opposite Sign Dilepton Events 

Cross section contours for the OS dilepton signal are shown in Fig. ^for the same cases as 
for the 1£ signal in Fig. ^. The solid lines are for = 100 GeV for which the SM background 
is 630 fb, while the dashed lines are for = 200 GeV for which the background is just 
9 fb. As in Fig. |^ and in subsequent figures, we show the dashed contours only for relatively 
large values of mi/2 for which employing the larger E^ value is really essential. A striking 
feature of Fig. ^ is the sharp kink near mo ~ 400 GeV where the contours change their 
slope. In cases (b) - (d) this is simply due to the opening up of the two body decays of 
the chargino and Z2 into u and iL (their branching fractions to are strongly suppressed 
because Z2 (Wi) has very small (zero) U{1) gaugino components). In case (a) however, the 
kink in the 1 and 2 fb contours occurs at around mo = 500 GeV. We have checked that 
this is because Z2 and Wi leptonic three-body decays mediated by left-handed sleptons have 
significant branching fractions (few percent) even though the two-body decays Z2 — >■ ZZi 
or Z2 H^Zi and Wi — > WZi are kinematically accessible: the resulting enhancement of 
the leptonic branching ratio, especially of Wi, accounts for the kink being somewhat beyond 
the sleptonic two-body decay region in case (a). 
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The 5cr observability level is at 40 fb for = 100 GeV and at 4.7 fb for = 200 GeV. 
We thus see from Fig. ^ that with = 100 GeV, the LHC should be able to observe a 
signal in this channel if mi/2 < 300 - 400 GeV (200-300 GeV, if we also require | > 0.25). 
The reach improves to mi/2 =400-500 GeV if the analysis is done using E'^ = 200 GeV. 
Notice that the reach is slightly larger in the tan /5 = 10 cases than in the low tan (3 cases (a) 
and {b). This is because the branching fraction for the two-body Z2 ZZ\ decay, which 
is very small for cases (a) and (6), is sizable when tan/? is large. We have checked that the 
statistical significance of the signal is marginally improved with E'^ = 300 GeV, but the 
cross section is then just around 1 fb for mi/2 = 500 GeV. In summary, with suitable cuts 
and 10 fb~^ of data, LHC experiments should be able to detect a signal in the OS dilepton 
channel for mi/2 up to 400-500 GeV, which corresponds to a gluino mass just beyond 1 TeV. 

If an OS dilepton signal is seen, one may again attempt to localize the position in 
parameter space via a measurement of the total OS dilepton cross section, which should place 
one along one of the contours in Fig. |^. Since the Z2 branching ratio into Z bosons depends 
on tan /5, the number of reconstructed Z ^ ii events may offer some rough discrimination 
in that parameter if the spoiler modes are kinematically accessible. An idea of mi/2 from 
the l^j, or li channels together with the cross section in this channel would enable the 
determination of mo if it is small: for example, for mi/2 = 400 GeV, the OS cross section 
rapidly varies from > 128 fb (small mo) down to 32 fb (mo ~ 500 GeV), and then slowly 
decreases to ^ 8 fb. Ratios such as a{OS)/a{li) or a{OS)/a{0i) would presumably be more 
accurately calculable than the absolute cross sections. 



It was noted in Ref. fT^ that the production of neutralinos in SUSY events can lead to 
a flavor asymmetry in the OS dilepton event sample, which may allow further parameter 
space location. For instance, if OS dileptons are primarily coming from Z2 decay, then they 
should mainly be of same flavor, e.g. ee or /i/i pairs. If instead OS dileptons come mainly 
from charginos or third generation quarks and squarks and their subsequent leptonic decays, 
then one would expect roughly equal abundance of e/2 and fie pairs as compared to same 
flavor lepton pairs. We have plotted the OS dilepton flavor asymmetry, 

A^(ee) + N{fifi) - N{efi) - N{fie) 



N{ee) + N{fiji) + N{epi) + A^(/ie)) 



in Fig. [ly for the same cases as in the previous figure. SM backgrounds are included in the 
figure. Points denoted by an x have an asymmetry Ap < 0.2, consistent with no asymmetry 
in our simulation. Open boxes or diamonds have asymmetry 0.2 < Ap < 0.5, while filled 
boxes or diamonds have Ap > 0.5. For mi/2 = 100 GeV, we took E^ = 100 GeV, and we use 
the box symbols. For larger mi/2, we took Efj^ = 200 GeV to improve the signal/background, 
and we use the diamond symbols. We see in frame (a) that there is a large asymmetry for 
mi/2 ^ 200 GeV, and also for small mo values. In the former case, for small mo, this is 
due to an enhanced branching fraction for Z2 iiZi, while in the latter case it is due 
in part to the Z2 £ii two-body decays. In frame (6), the asymmetry disappears for 
small mi/2 values due to interference effects driving the Z2 branching fraction to very small 



values pT| , p2| . The two frames for tan/3 = 10 continue to have significant flavor asymmetry 
even for mi/2 as high as ~ 300 GeV due to the significant Z2 ZZ\ branching fraction. 
Again the significant asymmetry for small values of mo is due to real slepton decays of the 
neutralino. Note also that in some cases there is an observable asymmetry even when the 
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leptonic branching fraction of Z2 is so small that the clean 3^ signal from W1Z2 production 
(discussed in Sec. IV) falls below the observable level. We thus see that an observation of a 
significant flavor asymmetry will localize us in the regions of the plane where at least one of 
mo or mi/2 is not too large. Furthermore, if event rates indicate a large value of mi/2, the 
observation of a flavor asymmetry would lead us to conclude that tuq is rather small. 

We also mention that we have checked that the jet multiplicity increases with mo for 
^1/2 = 200 and 400 GeV, with other parameters fixed as in Fig. |^. We have also checked the 
(ub) distributions which show a qualitatively similar trend as in the li case shown above. 
Again the results are essentially the same for the two signs of /i. We do not show these 
distributions for the sake of brevity. 



3. Same Sign Dilepton Events 



The SS dilepton plus jets channel has long been known to provide a clean 

signature for supersymmetry, and has been the subject of several studies P,|S|, |55|JT7| , |5D[] at 
the LHC Within the MSSM framework, detailed simulations have shown that, in this 
channel, the gluino reach extends to beyond 1 TeV. It is also known |17|J39(] that while gluino 
pair production with gluinos decaying via the chain g Wi — £ is frequently considered 
to be the main source of these events, many other sources may be important, in particular, 
decays of gluinos to third generation fermions and also squark decays. We also stress that 
same sign dilepton events do not necessarily originate via production of Majorana particles. 
For instance the production of bibi pairs, where 61 — > tWi can also lead to SS dilepton plus 
multijet topologies. It is clear that a reliable computation of the SS dilepton signal requires 
that all the decay chains as well as all possible production mechanisms be included, as is 
done in ISAJET. 

The SS dilepton cross section is shown in Fig. ^ for the same values of SUGRA param- 
eters as in Fig. H for = 100 GeV (solid contours) and E!j. = 200 GeV (dashed contours). 
As in Fig. ^ (and for essentially the same reasons) we see that the contours show a kink 
near the region where ii and u masses approach ~ "^Za" "^^^ backgrounds to the 
signal are just 1.7 fb and 0.25 fb, respectively, yielding "5cr" limits of 2.1 fb and 0.8 fb. (For 
E!j^ = 200 GeV, the Poisson probability of an expected background of 2.5 events fluctuating 
to 8 events is 4 x 10~^, so that the 10 event level is perhaps a more reasonable estimate of 



the reach.) We see from Fig. [IT] that even with E!j. = 100 GeV, the 5a reach extends out 
to mi/2 ~ 400 - 500 GeV, and up to 700 GeV in the small mo region, where squarks are 
relatively light and leptonic decays of Wi and Z2 are enhanced. The signal to background 
ratio exceeds unity. A higher value of E^ only gives a marginal increase in the reach. With 
the harder cut, the signal is small so that perhaps 20fb~^ of integrated luminosity may be 
necessary in this case. As in Fig. ^ the cross sections are somewhat larger in the tan /? = 10 
cases as compared to the tan/? = 2 cases. Finally, we note that although there are some 
fluctuations in our simulation for mi/2 < 200 GeV, this signal should again be observable 
down to relatively low values of mi/2- 

Again, a measurement of the total cross section for SS dilepton events will place us along 
one of the contours in the mo - mi/2 plane. As before, a measurement of the ratio of the 
SS cross section to the li cross section would be an indicator of the small mo, large mi/2 
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region. As in the OS dilepton case, the cross sections are somewhat larger for the high tan (5 
cases. It has been pointed out that a charge asymmetry may exist in any SS dilepton 
signal detected at a pp collider; this observation has since been confirmed by more detailed 
simulations 0. As for the asymmetry in the l^ channel, the SS dilepton charge asymmetry is 
again a reflection of the valence u and d quarks in the proton participating in the production 



mechanism. We show in Fig. [T^ the charge asymmetry, 

N{^+^+) - N{i-t 



N{£+i+) + N{i-t 



as a function of mo, for (a) mi/2 = 200 GeV (with = 100 GeV), and (b) mi/2 = 400 GeV 
(with = 200 GeV). We take Aq = 0, tan (3 = 2 and yU < 0. Frames (c) and (d) are same as 
(a) and (6), except for the sign of fi. It can be noted that for niQ large, where gg production 
is dominant, there is essentially no charge asymmetry. As tuq decreases, and squarks become 
lighter, the asymmetry grows, reflecting the presence of w-squarks as a significant source of 
the event sample. 

We have checked the dependence of the jet multiplicity in this sample on rriQ for the same 



cases as in Fig. |T2|. While we do see the anticipated trend for an increase of {rij) with mo, 
the extraction of mo appears more difficult than in the li case, in part because of somewhat 
larger error bars in our simulation. We have also checked the (ub) distributions for these 
same cases — they appear to be qualitatively similar to those for the li and OS dilepton 
samples. 

Finally, we have checked the sources of SS dilepton events for several points in parameter 
space to see how frequently these occur when each gluino decays via g — > qqWi (where q ^ t) 
and Wi decays leptonically since this chain has been suggested as a way for extracting 
rrig. For small values of mo — mi/2 this does not happen because gluinos decay to real 
squarks. For mi/2 = 250 GeV and mo = 4mi/2, we have checked the sources in cases (a) - 
(c) of Fig. We found that in these three cases, just 2/51, 3/36 and 0/48 events came 
from this particular cascade chain. Typically, in about half the SS dilepton event sample 
that passed the cuts, there was at least one t-quark from the decay of the gluino (recall 
that decays to third generation may be enhanced), and greater than a third of the events 
had their origin in gq^ or gqji production even though the squarks were somewhat heavier 
than gluinos. While these numbers do depend on the details of the cuts, they underscore 
the importance of simulating all possible production mechanisms and decay chains in order 
to make a realistic assessment of the feasibility of mass measurement in this channel. 



4-. Trilepton Events 

Finally, we show the cross section contours for 3£ + jets + ^j, events in Fig. again for 
the same four cases (a) - {d); as before, E^ = 100 GeV (solid) and 200 GeV (dashed). We 
note the following: 

• As expected, the cross sections are enhanced in the region where the two body decays 
of Wi and Z2 to and u are kinematically allowed. 

• The cross sections remain substantial even in the region of the plane where the spoiler 
decays of Z2 (the boundaries of these regions are denoted by dotted lines in Fig. ^ 
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become kinematically allowed. We thus conclude that while the cascade decay chains 
g —>■ qqZ2 — > qqilZi or g — > qZ2 qUZi are important sources of leptons in these 
events, there must be other sources operating as well. These include cascade decays 
of squarks and gluinos to charginos (which decay to W bosons) and to top quarks. 
We remind the reader that gluino decays to third generation quarks can frequently 
be enhanced [1^,^,^ because ti and hi are frequently lighter than other squarks (so 
there is more phase space), because the third generation Yukawa couplings can be large 
and enhance such decays, and because Higgs bosons that are produced in the decay 
cascades decay either to third generation fermions, or to charginos and/or neutralinos 
(which have enhanced decays to third generation fermions). 

• For mi/2 ^ 250 GeV (the boundary of the Z spoiler in Fig. |^) some of the trilepton 
sample (or the dilepton sample in Fig. |^) should consist of real Z + I (real Z) events. 

• There are regions of parameter space [0,^ where the leptonic decays of Z2, and 
hence the M signal, are suppressed. This causes a dip in the cross sections of frames 
(c) and {d) around (mg, mi/2) ~ (300, 200) GeV. It is instructive to compare this with 
the corresponding case for the clean trilepton signal from W1Z2 3i discussed in the 
next Section. 

We find background cross sections of 4 fb and 0.07 fb, respectively, for the two choices 
of E!j.. The 5a level for corresponds to a cross section of 3 fb, while the 5-10 event 
level might be a reasonable estimate for the reach with the larger value of i?^. We thus see 
that with Ej^ = 100 GeV, the reach in the trilepton channel extends up to 350-500 GeV 
depending on the parameters, except of course in the small mo region where mi/2 values as 
high as 700 GeV may be probed. Again, the reach is larger in the tan = 10 cases. With 
E^ = 200 GeV, the reach in mi/2 increases by about 50 GeV. However, since the background 
is essentially negligible in this case, it may be possible to push the limits even further with 
a larger data sample. 

As with the other multilepton channels discussed earlier, there should be distinctive 
properties of the 3i + jets + I^j, signal that allow some localization of where one is in 
parameter space. These include jet multiphcity, B multiplicity and various px distributions. 
Also, there should again be a charge asymmetry, where we expect more + H — events than 

H . We have also checked the variation of {rij) with mo. The distributions are somewhat 

flatter in this C3iSG clS compared to the li sample in Fig. |^ This is reasonable since for 
higher lepton multiplicity, the number of "partonic jets" must be correspondingly reduced. 
In contrast, the {ub) distributions are qualitatively similar to the previous cases. We do not 
present these plots here for brevity. 

5. A recapitulation of the LHC reach via multijet plus multilepton events 

We have seen that, regardless of the model parameters, for nig ^ 1 TeV (mi/2 ^ 400 GeV) 
there should be an observable SUSY signal in each of the li, OS, SS and 3i channels if 
the SUGRA framework that we have adopted is a reasonable description of nature. In our 
previous study [0], we saw that there will also be a clearly observable signal in the jets + ^^ 
channel. Hence, in SUGRA, a wide variety of supersymmetric signals are expected to occur 
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at the LHC If gluinos are heavier than 1 TeV, the signals in the di- and tri-lepton channels 
may not be observable, although signals in the multijet+^j, and 1£ channels may still be 
visible. The single lepton channel yields the maximal reach. Our computation shows that 
at the LHC experiments will probe mi/2 up to 600 - 700 GeV (m^ up to 1500 - 1800 GeV) 
even if mo is very large; if mo is relatively small, it will be possible to search for gluinos 



heavier than 2 TeV 42 



Before closing this discussion, we should also mention that we have examined the Ai 
channel 0. We find that, with = 100 GeV, these signals might be observable when 
"^■1/2 ~ 300 GeV for tan/? = 2 (or 500 GeV for tan/5 = 10), the reach in this channel is 
always smaller than in other channels. For this reason, and because there are rather few 
events in our simulation, we do not show these here. It should, however, be kept in mind 
that sparticle production can lead to these striking events at an observable level, especially if 
mo is small and mi/2 not very large. Reducing or even eliminating the jet cuts could lead to 
larger signals in these event topologies without large increase in the background (assuming 
that leptonically decaying Z bosons can be readily identified). We do not consider this any 
further in this study. 

Up to now, we have fixed Ao = in our analysis. The cross sections should mainly 
depend on Aq due to the variation of third generation squark masses. Instead of performing 
lengthy scans of the parameter space, we have illustrated the Aq dependence of the cross 
section in Fig. for six choices of Aq and for mo = 500 GeV, mi/2 = 160 GeV, tan (3 = 2 and 
/i > 0. For y4o = 0, 500 and 1000 GeV, gluinos decay via three-body modes into quarks plus 
various charginos and neutralinos. For larger, positive values of Aq, the decay patterns of ti 
are qualitatively similar to those for Aq = 1000 GeV until Aq exceeds ~ 1330 GeV, at which 
point m?^ becomes negative. For the three negative values of Aq sampled, the top-squark is 
so light that g —>■ tti dominates the gluino decay channels (this is sensitive to other model 
parameters, including sgn yu). For Aq = —500 and —900 GeV, the ti dominantly decays via 
ti — > bWi, while for Aq = —950 GeV, it is so light that only ti —>■ cZi is allowed. For all 
Aq values sampled, we see that the li cross section is roughly constant to within a factor of 
< 2. The dilepton cross sections show a somewhat larger variation, although this may not be 
sufficient, by itself, to determine Aq, since small changes in mi/2 can cause similar variation. 
The maximum variation is seen in the 3i cross section. The dilepton and trilepton cross 
sections are largest for cases with large negative Aq values of for which g tti and g tti 
are the only two body decays of the gluino. Since ti decays via ti — > bWi, the increase in the 
leptonic cross sections (particularly for the 3i channel) should not be surprising. The sharp 
drop in the cross sections at the most negative value of Aq is because the chargino decay 
mode of the t-squark (which is a source of leptons) becomes inaccessible, and ti cZi. In 
this case, gg pairs (with g tti) can give rise to events with at most two hard, isolated 
leptons, and a{SS) ~ a{OS). It may ultimately be possible from the ratios of multilepton 
to single lepton cross sections to pin down Aq, especially if it is close to the boundary of the 
excluded region where becomes negative. We have, however, seen that this ratio shows 

a similar trend in the small mo region where Wi and Z2 leptonic decays are enhanced. A 
measurement of (ub) could serve to distinguish the two different origins of leptonic signals. 
We also note that in principle, there could be parameter values for which g bbi might be 
the only allowed two body gluino decay, in which case we would expect a reduction of the 
multilepton cross sections. More detailed study of the variation of the signals with Aq are 
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clearly necessary before definitive conclusions can be drawn. 



Our preliminary conclusions based on Fig. are that i) the multi-lepton cross sections 



examined above are less sensitive to variation in Aq for lower lepton multiplicity, and ii) 
except for the extreme cases where new channels for gluino decays open up (these might be 
signaled by events with unusually high i?-hadron multiplicity), even the multilepton cross 
sections are rather insensitive to Aq, and the choice Aq = that we have adopted yields 
representative values of these cross sections. 

IV. CLEAN MULTILEPTON SIGNATURES FOR SUPERSYMMETRY 

In the previous section, we focussed on the study of multilepton events with at least 
two hard jets and substantial The cascade decays of gluinos and squarks were the 

main source of these jetty events. While the direct production of charginos, neutralinos and 
sleptons can also lead to similar event topologies, these signals would be more difficult to 
pick out from SM backgrounds because of relatively lower total cross sections and softer 
PrU^t) and I^j, distributions. Moreover, there would be the additional issue of how to 
separate them from the corresponding signals from gluino and squark cascades for which 
the cross sections are considerably larger. Clean multilepton events, i.e. events without any 
jet activity, for which SM backgrounds are smaller, offer a more promising way of searching 



for chargino and neutralino |26,Tq] or slepton |24,25| signals at the LHC We study the reach 



in SUGRA parameter space in these channels in this section. 

Unlike in the previous section, where for each point in SUGRA parameter space we gen- 
erated all SUSY subprocesses using ISAJET, here we focus on specific sets of reactions. This 
is because the majority of events generated contain gluinos and squarks which almost always 
yield hard jets, so that the efficiency for generating clean multilepton events is very small: 
the computer time that would be necessary to obtain an adequate sample of clean multi- 
lepton events would then make global scans of SUGRA space quite intractable. For sample 
points in the parameter space, we have checked how various SUSY channels contribute to 
the specific reactions that we are searching for. 

A. Clean Trilepton Events from W1Z2 Production 
These signals have previously been studied P^,[TB| within the framework of the MSSM 



for parameter sets motivated by SUGRA models. In our previous study [jl6|, we had fixed 
/i = —rrig and chosen tan (3 = 2 and = + GeV. We found that it was possible to find 
cuts which not only reduce SM backgrounds to negligible levels, but also isolate trileptons 
produced via pp — >• W1Z2 + X — (iuZi + I'I'Zi + X from those produced by other SUSY 
reactions. These other SUSY processes typically contribute ^ 10% of the total trilepton 
signal, at least for the parameters where the signal was deemed to be observable. Here, 
we extend our previous study and explore the reach of the LHC for this signal within the 
SUGRA framework, and delineate the region of parameter space where the clean trilepton 
signal should be observable above SM backgrounds. 
Exactly as in Ref. |T^, we require: 



three isolated leptons, with Pt(^i,4) > 20 GeV, Pt(4) > 10 GeV. 
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• a) a central jet veto, i.e. no jet with PT(jet) > 25 GeV within \rij\ < 3; 

• zzz) < 100 GeV; 

• iv) \m{ii) — Mz\ > 8 GeV for all combinations of OS leptons with the same flavor in 
the trilepton event. 

Cuts ii) and Hi) greatly reduce the backgrounds from the cascade decays of gluinos and 
squarks, while iv) is designed to eliminate WZ events. After these cuts tt remains the 
dominant background. It can be greatly reduced by further requiring, 

• v) the two fastest leptons have the same sign of charge and the flavor of the slow lepton 
be the anti-flavor of either of the two fast leptons. 

This reduces the signal by 50% but essentially eliminates the top background, from which the 
two hardest leptons almost always come from the primary decays of the t-quarks, and hence, 
have opposite signs of charge. To recover some of the rejected signal without a significant 
increase in the ti background, 

• vi) we retain events in which the two fastest leptons have opposite sign provided 
prih) > 20 GeV. 

After cuts i) - iv) and either v) or vi), we find a SM background level |[10| of 0.7 fb from 
WZ production where the gauge bosons decay into e, fioi r (which then decays leptonically), 
and 0.13 fb from ti production (for rrit = 170 GeV), yielding a total SM background of 



0.83 fb. This is somewhat larger than in our earlier study because of differences in parton 
distributions as well as calorimeter simulation. Assuming an integrated luminosity of 10 fb~^. 



the minimum signal cross section for observability at the "5(T level" (Ai'signai > 5yA'bkgds) 
works out to be 1.44 fb: the Poisson probability of an upward fluctuation of this amount is 
2x10"^ Notice that iV^ignai/iVbkgd > 1-7. 

The region of the mo - mi/2 plane where the signal is observable at the 5a (lOa) level is 
shown by hollow (solid) squares in Fig. ^ for Aq = and (a) tan/3 = 2, /i < 0, (b) tan/3 = 
2,/i > 0, (c) tan/? = 10, /i < 0, and (d) tan/3 = 10, /i > 0. For each parameter space point 
sampled, we require at least 25 events to pass the cuts in our simulation. The x's show the 
points that we have sampled but for which the signal falls below the 5a level. Also shown in 



Fig. |15| are the boundaries of the region where the spoiler modes Z2 — > ZiH^ or Z2 —>■ Z\Z, 
or two body lepton-slepton decays of the neutralino become accessible. In cases &) - d) the 
boundary of the Higgs spoiler decay is not shown as it always lies above the boundary of 
the Z spoiler. Several features of this figure are worthy of mention: 



In case a), which corresponds most closely to the points sampled in Ref. [|T6|], we see 
that the signal is observable at the 5a level all the way up to the boundary of the 
spoiler modes, and for most of the region the significance is larger than lOa. 

There are regions of the mo - m.1/2 plane in cases h) - d) where the chargino is at its 
current experimental bound from LEP, but where the trilepton signal fails to satisfy 
our 5a criterion for observability. This was traced [p6i^p2| directly to the leptonic 
branching fraction of Z2 which can drop by as much as two orders of magnitude due 
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to interference effects between the slepton and Z mediated decay amplitudes. Thus a 
non-observation of a signal in this channel will not allow us to infer a lower limit on 
either or vn^^. The regions of the parameter plane where there is an observable 
signal in this channel at the LHC are similar to the regions that the Tevatron operating 
at 10^'^cm^^s^^ could probe 



Except in the "hole" mentioned above where there is no observable signal, the trilepton 
signal should be detectable all the way up to the limit of the spoilers. If sleptons are 
light enough so that Z2 — > ^l,r^ are kinematically accessible (the small mo region of 
the plane), then these decays may dominate the spoiler decays. Then the branching 
fraction for leptonic decays of Z2 is very large, and the reach in the trilepton channel 
extends well beyond the boundary where the spoilers become accessible. Notice the 
small wedge between the contours labeled Z2 — >■ and Z2 — *■ vv where the signal 
drops because the invisible decay Z2 — > vv of the neutralino dominates. 

We see that flipping the sign of /i makes a much larger difference in the tan 13 = 2 cases 
a) and h) relative to the tan (3 = 10 cases c) and d). This can be understood if we 
recall that it is always possible to choose /i and the gaugino masses to be positive by 
convention: then, the vacuum expectation values of the two Higgs fields can no longer 
be chosen to be always positive, and the physically relevant sign between fi and the 
gaugino masses appears as the sign of tan /?. Of course, for large values of tan P where 
one of the vacuum expectation values is essentially negligible, this sign is unimportant, 
explaining why the results in cases c) and d) appear so similar. 



MSSM case studies of Ref. [|T6| suggest that the clean trilepton signal is relatively pure, 
and that the "contamination" from SUSY sources other than W1Z2 production is small. It 
should be kept in mind that in these studies we had fixed tan/3 = 2 and chosen fi = —rrig, 



so that the situation is roughly that in Fig. Oa, where the signal exceeds lOa over most of 



the plane. There are substantial regions of the parameter plane in cases b) - d) where the 
significance of the signal is between 5 and lOa. 



It has already been pointed out [|T6| that the isolation of the signal from W1Z2 production 
will allow a reliable determination of m^^ — m^^, and perhaps also other combinations of 
chargino and neutralino masses. 



B. Clean Dilepton Events from WiWi Production 

We have just seen that while charginos and neutralinos might be detectable over a large 
regions of parameter space in the clean trilepton channel, there are parameter ranges for 
which the leptonic decays of Z2, and hence this signal, is strongly suppressed even if charginos 
are relatively light. We are thus led to examine whether OS dilepton signals from the reaction 
pp — > WiWi + X — > IvZi + ^'vZl + X might be able to probe charginos in these regions, 
or to provide a new channel for confirmation of the existence of charginos detected in the 
trilepton channel; this was found to (at least partially) be the case for the 10^^ cm~^s~^ 



upgrade of the Tevatron |^ 



To search for events in this channel we have made the following cuts: 
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• We focus on e^/i^ events with |pt(^)| > 20 GeV to eliminate large backgrounds from 
Drell-Yan production; 

• We veto events with any jet with Et > 25 GeV within |?7| < 3; 

• We require 30° < A0e^ < 150°; 

• We require A0(pr(e/i), -^7.) > 160°; 

• We require 40 GeV < J^rp < 100 GeV (the upper limit on I^j, is to prevent other SUSY 
sources from contaminating the signal). 

We have used ISAJET to compute SM backgrounds to the dilepton signal from tt, WW, 
rf, WZ and ZZ production. We find that our cuts efficiently suppress backgrounds from 
all but WW events, for which the cut cross section is 136 fb (compared to the a(ti) = 9.9 fb 
and (j(rr) = Ifb). For an integrated luminosity of 10fb~^, the 5a level of observability 
corresponds to a signal cross section of 19 fb, although the signal/background ratio is small. 
We sampled points in the mo - mi/2 plane for the same cases as for the clean trilepton signal 
in the previous subsection. We found that except for a few points near mo = in case a) 
and an isolated point in case c), the signal is below the 5a level, and for most of the plane, 
even below the 3a level. We conclude that, unlike at Tevatron upgrades, the dilepton signal 
from chargino pair production is unlikely to be observable above SM backgrounds. 



C. Clean Dilepton Signals from Slepton Pair Production 

Charged sleptons and sneutrinos can be pair produced at the LHC in qq fusion processes 
via charged or neutral gauge boson exchange in the s-channel. Their (cascade) decays can 
lead to event topologies with several leptons and jets in the final state. Previous studies 
2^,25] have shown that the clean, acollinear e~^e~ + and fi~^fj^~ + I^rp channels offer the 



best hopes for the discovery of sleptons at the LHC. Our main purpose here is to delineate 
the region of the SUGRA parameter space where these signals might be observable at the 
LHC, and to check whether these can be distinguished from corresponding signals from 
chargino pair production. 

To separate the signal from SM backgrounds, we require ||25|| , 



• i) exactly two isolated same flavor OS leptons, each with |pt(^)| > 20 GeV, 

• 11) > 100 GeV, 

• in) a veto on central jets with Et > 25 GeV within |?7| < 3, and 

• iv) A(p{pT{ii),fT) > 160°. 

After cuts i) - iv), the dominant SM backgrounds to the SUSY signal come from ti 
(2.2 fb) and W+W- (2.9 fb) yielding a "5ct observability level" of 3.6 fb for a year of LHC 
operation at the design luminosity. The slepton cross section is, however, rather small and 
a higher reach is obtained with somewhat stiffer cuts to further reduce the background at 
modest cost to the signal. Hence, 
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• v) for detection of heavy sleptons, we also require |pt(^)| > Pt A(f){ii) < A0c, 
where and A0c can be adjusted appropriately. In our analysis, we fix = 40 GeV 
and A0, = 90°. 

Including cut v), we find no events pass the cuts from our simulation of the WW sample 
(a one event level corresponds to a = 0.0015 fb in our simulation) while from tt events, we 
find a background cross section of 0.07 ± 0.006 fb. We thus expect < 1 background event 
per LHC year, with our "hard" cuts i) - v). 

The region of the mo - rni/2 plane where slepton production should yield observable 
signals after these hard cuts is shown in Fig. |l^, again for the same four cases as in previous 
figures. Since the SM background level is very small, we show contours of constant cross 
sections corresponding to the 3-5 event level per LHC year (triangles), 5-10 event level 
(hollow squares), 10 - 20 event level (squares with crosses) and > 20 event level (filled 
squares). The crosses denote the sampled points for which the cross section is smaller than 
3 fb. We also show contours where m^^ = 70, 100, 150, 200 and 250 GeV. If we take the 
five event level to give the optimistic reach, we see that the reach of the LHC extends to 
~ 250 GeV, corresponding to m^^ and m,> to just over 300 GeV for larger values of mi/2- 
For a SM background expectation of 0.7 events, the Poisson probability of a fluctuation to 
the five (ten) event level is 8 x 10^^ (4 x 10^^), so that a conservative estimate of the reach 
after a year of LHC operation is somewhere between 5 and 10 events. 

More disturbing is the existence of the "hole" where the cross section falls below the five 
event level for small values of mo and mi/2 in cases a) and c). Notice that unless the energy 
of LEP2 is upgraded so as to ensure the detectability of sleptons as heavy as 100 GeV, 
ifi (and, of course, also u and £l) may evade detection at both LEP2 as well as at the 
LHC. To understand why the hole is much larger for the /i < cases, we have examined 
the differences in sparticle properties for {niQ, 1711/2) = (40 GeV, 140 GeV) in cases a) and 
b). For the negative case a), Wi and Z12 are somewhat heavier than in case b) so that 
the mass difference between and Zi is rather small (14 GeV, in our example). As a 
result, the efficiency for particularly inin events to pass the hard Pt{^) > Pt and cuts 
is reduced, leading to a drop in the cross section. For case b) nii^ — m^^ = 30 GeV so that 

the daughter leptons are considerably harder. In addition, II predominantly decays to Wi 
and Z2, and further, the leptonic branching fraction for Wi is enhanced to 22%, while the 
neutralino decays via Z2 ini, so that hard leptons can come via several chains. This 
example also underscores the importance of incorporating the various cascade decays into 
the slepton analysis. 



We have just seen that because of the hard cuts that we have used in Fig. |16], there are 
small regions of parameter space where sleptons with masses ~ 80 - 120 GeV may evade 
detection both at LEP2 and at the LHC. Because of the importance of this issue, we have 
re-done our analysis using just cuts i) - iv) for which the SM background cross section is 



5.1 fb. In Fig. we show the regions of the mo-mi/2 plane where the significance of the 



signal a = Ai'signai/-\/-^bkgd is 3(j (triangles), 5a (hollow squares), lOcr (squares with crosses) 
and 20(7 (filled squares) for the same four cases as in Fig. |16|. Indeed we see that with the 
soft cuts, the slepton signal always exceeds the 5a limit in the "hole" regions of Fig. |T^, 
and further, that there is no window of masses where sleptons will escape detection both 
at LEP2 and at the LHC. The maximal reach at the LHC is, of course, obtained using the 
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hard cuts. 

In order to check whether dilepton events from slepton pair production might be confused 
with corresponding events from chargino production, we have checked the origin of the events 
which satisfy our cuts for several cases: 

• hard cuts, with (mo,mi/2)= (210 GeV, 160 GeV), case a) for which the slepton masses 
are 225 - 250 GeV and m^^ = 155 GeV; 

• hard cuts with {mo,mi/2)= (60 GeV, 160 GeV), case b) for which the slepton masses 
are 92 - 135 GeV and m,j, = 109 GeV; 

Wi ' 

• hard cuts with (mo,mi/2)= (40 GeV, 140 GeV), case b) for which the slepton masses 
are 75 - 115 GeV and m,j, = 86 GeV. 

Wi 

In all these cases, although we had generated all slepton (including sneutrino) as well 

as WiWi events using ISA JET, we found that only slepton events in our sample of 40 - 
60 events that pass our cuts i.e. there were no events from direct chargino pair production 
in the sample. We did find events from cascade decays of sneutrino (produced in pairs 
or along with a charged slepton) as well as fs. To check whether chargino production 



contaminates the slepton sample with the soft cuts i) - iv) in Fig. [T^, we have checked the 
sources for case a) with (mo,wii/2) = (40 GeV, 140 GeV), for which the cross section is 
< 0.3 fb after hard cuts, but where the signal exceeds 5a with soft cuts. We find that out of 
a total of about forty events that pass the cuts in our simulation, just six come from direct 
chargino production, with the charginos decaying via iui, {i = e,/i). We thus conclude 
that a conclusive observation of a dilepton signal with the hard cuts will be unlikely to 
be confused with chargino pair production. There may, however, be some small chargino 
contamination of the signal with the soft cuts. In this case, the event sample should be 
large enough to provide other handles on chargino-slepton discrimination. For instance, if 
rriej^ = rrijy,^ chargino production should lead to as many e^^"^ events as e^e~ events, 
whereas we would expect significantly more same flavor events in the case of slepton pair 
production. 

Finally, for the first two cases with the hard cuts above, as well as for the the soft cut case 
we just discussed, we generated all SUSY subprocesses and ran them through the "slepton 
cuts" to see whether the "slepton signal" is contaminated by squark and gluino production, 
which occurs with much larger cross section. This requires a simulation of a very large 
number of events since only a very tiny fraction of events pass the cuts. We examined the 
twenty events that satisfied the "slepton cuts" in each of these three cases: we found just one 
event from squark and gluino sources in one of the three event samples. However, in almost 
half the events for the /i < 0, hard cut and the soft cut cases, the leptons both originated 
from Z2 decays in W1Z2 or Z2Z2 events (the leptons from the decay of a single Z2 satisfy the 
A(f) cut more readily than those from H^iW^i events). In the yU > 0, hard cuts case simulated, 
we have m^^ > rn^^^^, m^^ > fni^, so that Z2 always decays via Z2 Hr into real sleptons: 

for this case, we found about 80% of the events had their origin in Wi and Z2 production. 
We thus conclude that while squark and gluino production is unlikely to contaminate the 
slepton sample, Z2 decays from W1Z2 or ^2^2 production can significantly contaminate the 
slepton signal (presumably WiWi events frequently fail the Acj) cut, which fails to remove 
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dileptons from Z2 decays). However, these processes will themselves lead to characteristic 
signatures (the clean trilepton signature discussed above or even 4£ topologies) and would 
be detectable in their own right. We should also add that since these have not been included 
in Fig. 16 and Fig. 17, the actual cross sections may be somewhat larger than shown in these 
figures. 

Before drawing final conclusions regarding the detectability of sleptons at the LHC, we 
stress that we have assumed a 100% jet rejection efficiency for jets in the fiducial region. A 
real detector will, of course, have cracks and other dead regions. This is especially important 
here because the crucial cut for the detectability of sleptons over the background from 
tt production is the central jet veto. In our previous analysis p5|, we had shown that with 
the hard cuts, the tt background increases by about a factor of about five if instead this 
veto efficiency is 99%. Except to point out that it may be possible to reduce this detector- 
dependent background significantly by adjusting and A(/)c, we will not discuss this any 
further. We thus conclude that if detectors have the capability to veto central jets with a 
high efficiency, it should be possible to probe l^i and /Ir masses up to about 250 GeV at the 
LHC. The slepton signals are, however, very small so that perhaps 20 - 30 fb^^ of integrated 
luminosity may be necessary to confidently probe their existence. 

Finally, we point out that in a recent paper [Q, the cosmological relic density from 



neutralinos produced in the early universe was evaluated for the same SUGRA model. In 
these calculations, it was found that a relic density of ^Ir? ~ 0.15 — 0.4, which is favored by 
cosmological models with a critical density and a 2:1 mixture of cold/hot dark matter, would 
occur mainly if the slepton mass m^-^ ~ 100 — 250 GeV. Thus, failure to detect a slepton at 
LHC could place rather severe constraints on cosmological scenarios which ascribe the bulk 
of cold dark matter in the universe to stable neutralinos. 



V. COMPARISON OF RESULTS FROM THE VARIOUS CHANNELS 

We have used ISAJET to map out the region of parameter space of the minimal SUGRA 
model with radiative breaking of electroweak symmetry where various n-jets plus m-leptons 
(n > 2, m = 1, 2, 3) plus signals are observable above SM backgrounds at the LHC. These 
signals are dominantly expected to come mainly from gluinos and squark production followed 
by cascade decays. This paper is a continuation of our previous study |ll[] where we had 



focussed on multijet plus events with an isolated lepton veto to reduce backgrounds from 
vector boson and top quark production. We also examined the reach in the complementary 
clean dilepton and trilepton channels to investigate the detectability of the electroweak 
production of sleptons and charginos/neutralinos at the LHC. 

Since the parameter space of the model is rather large, it is impractical to sample all 
regions of this space. One approach would be to generate random sets of model parame- 
ters (mo, ""7.1/2, tan/5, Aq, sgn/x) and investigate various signals for the set of models thus 



obtained. This is the strategy used in Ref. ||22| where the authors generated about IK pa- 



rameter sets in their exploration of the SUSY reach of the Tevatron and its possible upgrade 
options. While this is indeed a viable strategy and may indeed have the advantage that it 
samples the parameter space "more uniformly", it has some shortcomings. First, one has to 
choose how to sample each direction; e.g. should one randomly generate mo or In mo, since 
the measure on parameter space is unknown. This is important because (for each sign of /i) 
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just IOOO4 = 5.6 points are generated on average along each of the four directions. Second, 
and more importantly in our view, while it is true that there may well be a fairer sampling 
of parameter space with this approach, it is difficult to relate the results to the underlying 
parameters of the theory. For these reasons, we have chosen to perform detailed scans in the 
mo - mi/2 plane (sparticle masses which dominantly determine the rates and distributions 
of the various signals are most sensitive to these parameters) for fixed values of tan/3 and 
Aq. We illustrate the results for a small (tan/? = 2) and a medium (tan/5 = 10) value 
of tan f3. We do not consider larger values of tan f3 because the effects of bottom and tau 
Yukawa interactions, which could become important, have not yet been completely included 
in ISAJET. In most of our analysis, we fix = (this does not mean that the weak scale 
value of the A-parameter vanishes) since our signals are moderately insensitive to this choice 



(see Fig. [IJ) except very close to the boundaries of the parameter space region where the 
correct pattern of electroweak symmetry breaking is not obtained. 

The details of our calculation in the multijet channels may be found in Sec. Ill, while the 
clean multilepton signals are discussed in Sec. IV. Instead of repeating this discussion one 
more time, we have chosen to summarize the results for the LHC reach in the various channels 



in Fig. pifor Aq = and (a) tan/3 = 2,/i < 0, (b) tan/3 = 2,/i > 0, (c) tan/3 = 10,/i < 0, and 



(d) tan /? = 10, /i > 0. As before, the hatched (bricked) regions are excluded by experimental 
(theoretical) constraints. For a signal to be regarded as observable 03, we require that for 



an integrated luminosity of 10 fb at the LHC, the event rates and numbers satisfy: 



• a statistical significance > 5cr, where a = Nsignsi/ ^jN^^kgd', 

• A'signal/A^bkgd > 0.2; 

• A^signal > 5. 

In the region Fig. |18| below the dashed line (labeled I^j.) the lepton plus I^j, signal should 
be observable beyond the 5a level for an appropriate choice of the cut variable defined 



in Sec. Ill as well as in Paper I [|IT| from which these contours have been taken. The various 
dashed-dotted contours mark the boundaries of the region where the 1£, same-sign (SS) 
dilepton, opposite-sign (OS) dilepton and trilepton (3£) plus multijet plus signals should 
be observable at the LHC, again for some value of < 200 GeV (400 GeV in the case of 
the li signal) as obtained from the analysis in Sec. HI. The regions below the dotted line 
(labeled i) and solid line (labeled W1Z2) are where the clean dilepton and trilepton signals 
are observable as discussed in Sec. IV. 
Several comments are worth noting: 

• At the LHC, it should be possible to detect gluinos as heavy as 1.5 - 1.8 TeV {nig ~ 
2.3 TeV if rriq ~ f^g), corresponding to mi/2 < 600 - 700 GeV, after just one year 
of running at its lower design luminosity option of 10fb~^/year. This is considerably 
beyond the bounds (mi/2 ^ 400 GeV) obtained from (admittedly subjective) fine- 
tuning arguments, and so should provide some safety margin for the detectability of 
SUSY at the LHC, at least within this minimal framework with conserved i?-parity. 
We also remark that we found no holes where these signals (or the multilepton signals, 
for that matter) might escape detection. 
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It is interesting that the maximal reach is obtained in the li channeL This is because 
there are numerous sources of leptons in SUSY events so that a lepton veto significantly 
reduces the signal cross section. Our analysis using the parameter shows that 
backgrounds from W boson and tt production (which lead to isolated leptons in the 
final states) can be controlled without vetoing events with leptons. It should thus be 
possible to combine the signals in the I^j, and li channels to obtain a somewhat larger 
reach. 

If squarks and gluinos are lighter than 1 TeV, several other signals should be observ- 
able above SM backgrounds if a signal in the I^j, or li channels is to be attributed to 
sparticle production. Although our conclusion, strictly speaking, has been obtained in 
the rather constrained SUGRA framework, including constraints from radiative elec- 
troweak symmetry breaking (these essentially fix previous analyses suggest 
that this will be true even if constraints from electroweak symmetry breaking are 
relaxed [Q. We also note that a portion of the multileptonic signals arise from lep- 
tonically decaying Z bosons. This is the reason why the reach in multilepton channels 



is slightly larger in the tan/3 = 10 cases (c) and (d) in Fig. |T8[ The real Z boson sig- 
nals are sensitive to the value of /x 0, and hence, to the radiative symmetry breaking 
constraint. 

While it appears that only a rather small subset of the parameter plane can be probed 
via the clean leptonic channels, the observation of these signals is important because 
it leads to direct detection of Wi, Z2 (this sparticle may be hard to detect even at the 



NLC) and the sleptons. Moreover, it has been shown [jT6[ that it is possible to isolate 
W1Z2 — >■ 3i events from SM backgrounds as well as from other SUSY sources. This 
allows for a reliable determination of m^^ — m^^, and perhaps, other combinations 
of chargino and neutralino masses. We stress that the non-observation of a trilepton 
signal at the LHC will not lead to a bound on the chargino or neutralino mass because 
of parameter space regions where the leptonic decays of Z2 are strongly suppressed. 
It is, however, interesting to note that even in these regions, the multi-jet plus 3i 
signals are observable, implying that there are significant other sources of leptonic 
events (notably, third generation fermions and sfermions). 

At the LHC, it should be possible to detect sleptons with masses up to 250 GeV 
(300 GeV for ii) in the clean OS dilepton channel. We have also shown that sleptons 
as light as 80 GeV ought to be detectable at the LHC using the "soft cuts" discussed 
in Sec. IV. Thus, there is no window where sleptons might escape detection, both at 
LEP2 and at the LHC. Furthermore, the LHC is sensitive to the most favored range of 
slepton masses expected from calculations of the dark matter neutralino relic density 
(m^- ~ 100 - 250 GeV) 



Aside from the question of the detection of SUSY, it is interesting to ask whether it is 
possible to devise tests of the various assumptions underlying the minimal SUGRA frame- 
work that we have adopted for our analysis. Tests that work well at an electron-positron 
collider |^,^ do not appear to be feasible at the LHC, partly because the initial state of the 
colliding partons is not known, and partly because of the messy interaction environment at 
the LHC. Alternatively, we may ask whether it is possible to use the multitude of observables 
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that should be accessible at the LHC to determine the underlying parameters of the model. 
This is clearly a complex task since the directly observable quantities such as cross sections 
in various channels depend on various masses and mixing angles which have to be unraveled 
in order to get at the underlying parameters. 

In this paper, we have made a first attempt to understand how it might be possible 
to use the LHC data to get at mo and mi/2- We have little to say at present about the 
determination of tan (3, Aq or sgn /x. 

1. If mi/2 ^ 300 GeV (so that gluinos are lighter than about 700 - 800 GeV) we had 
shown in Paper I that it should be possible to measure rrig to 15-25% by requiring 
hemispheric separation of events in the channel. Presumably, the same strategy 
can also be used in the li channel. The value of rrig can be directly related to mi/2 
(aside from the (usually small) corrections due to differences between the running and 
pole gluino masses). 

2. If the trilepton signal from W1Z2 production is observed at a substantial rate, it 
would be possible to check whether the value of m^^ — m^^ is in agreement with 
the expectation from the gluino mass, assuming that is large and the unification 
condition for gaugino masses is valid. If the gluino and neutralino masses are not 
in accord with this expectation, we would probably conclude that |/i| is not large 
which would imply that we are somewhat close in parameter space the boundary of 
the bricked region where the correct pattern of electroweak symmetry breaking is not 
obtained. The alternative would be that the gaugino mass unification condition is 
invalid. 

3. If the various multijet signals are observed at rates compatible with gluino masses 
corresponding to mi/2 < 200 GeV, but no clean trilepton signal is seen, we would 
probably infer that we are in one of the "hole" regions where the leptonic decays of Z2 
are strongly suppressed. This would imply that mo, and hence, squarks and sleptons 
cannot be too heavy (although there would be no guarantee that sleptons would be 
light enough to be observable.) 

4. If a signal is observed in the OS dilepton channel with the "slepton cuts" of Sec. IV, we 
would place ourselves in the bottom left corner below the dotted line in the mo-mi/2 
plane. In this case the multijet topologies from gluino and squark production must be 
seen. Otherwise, the assumptions of universal sfermion and/or gaugino mass at the 
ultra-high scale, which imply m| = m| -|- (0.7 — 0.8)m? could not be valid. 

It is possible that gluinos are rather heavy so that neither the W1Z2 nor the slepton 
signals are accessible. The determination of parameters is more difficult in this case. We 
have shown in Sec. Ill that the cross sections for multijet plus lepton signals will place us 
on one of the contours in Fig. |^, |^, ^ or Because the multilepton contours are roughly 



horizontal (except in the rriQ < 400 - 500 GeV region, to which we will come back to), it 
should be possible to get a rough idea of mi/2 (roughly within ±(50 - 100) GeV) and hence, 
of rrig. It should also be possible to decide whether mo is small (^ 300 — 400 GeV) or 
rather large, with a degree of confidence by studying the ratio of a{Oi + jets) /a{ni + jets) 
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for n = 1 - 3. For small values of mo and somewhat large values of mi/2, the leptonic 
decays of charginos and neutralinos, and hence the multilepton signals, are enhanced. For 
the same reason, the Oi + jets signal (because of the lepton veto) is reduced, as can be seen 
by the down turn of the corresponding contour in Fig. |18|. This could be confirmed by a 
measurement of the flavor asymmetry in the OS dilepton sample (see Fig. |I0]). If gluinos are 
heavy, and mo > 500 GeV, the determination of mo may be more difficult. Possible handles 
are the charge asymmetry in the li and SS event samples (the asymmetry reduces with mo) 
or the jet and B multiplicities in the 0, 1 and 2 lepton multijet samples (the multiplicity is 
larger for larger values of mo) . Clearly detailed case studies beyond the scope of this analysis 
would be required to determine how well these model parameters can be determined. 

We have not found any strategies for the determination of tan /5, Aq or sgn /i. A qualita- 
tive idea of whether tan f3 is small (close to unity) or large might be obtained by looking for 
multijet events with real Z bosons: these are more abundant for larger values of tan f3. The 
observation of the Higgs boson and a measurement of its mass (perhaps in the 77 channel) 
may also provide a handle on this parameter: since rriHi = at tree level if tan/3 = 1, the 
lightest Higgs boson tends to be lighter when tan/5 is close to unity. The parameter ^0 
mainly affects the third generation. Variations in Aq can alter significantly the dominant 
gluino decay channels, so that rates in di- and, especially, trilepton plus multijet channels 
can have significant dependence on this parameter. The multiplicity of central 5-hadrons 
in SUSY events should also be sensitive to Aq. We have, however, not studied this aspect 
of parameter space in enough detail to draw any clear conclusions, on Aq. 

It may well be that all the parameters will ultimately be extracted by a global fit to all 
the data. The success of such a fit would certainly be non-trivial since the complete set of 
observations would need to be fitted by just four parameters (plus a sign). If an adequate 
fit is not possible, the assumptions underlying the model would need re-examination. 

To sum up, if supersymmetry is the new physics that ameliorates the fine tuning problem 
of the SM, it appears almost certain that there will be a multitude of new physics signals at 
the LHC Although our analysis has been performed within the framework of the -R-parity 
conserving minimal SUGRA model, we do not expect the results to be qualitatively altered 
due to minor modifications of the model, as long as i?-parity is conserved. The maximal 
reach is obtained in the single lepton channel and it appears that gluinos as heavy as 1.5 
- 1.8 TeV (2.3 TeV if squarks are degenerate with gluinos) ought to be detectable at the 
LHC with just 10fb~^ of data. It should also be possible, in at least some cases, to identify 
the sparticle origins of various signals. We have also made a preliminary exploration to see 
how one might attempt to localize the underlying SUGRA model parameters, given that 
these SUSY signals are seen at the LHC. While this may well be easier at e~^e~ colliders 
(with sufficient center of mass energy), it is certainly worthwhile to think about what might 
be possible in experiments at the LHC, where construction has already been approved. We 
have argued that it might be possible to extract mi/2 and, to some extent, also mo via 
a simultaneous study of several signals. Other parameters appear even more difficult to 
obtain, but this study should only be regarded as a first attempt in this direction. 
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TABLES 

TABLE I. Results of background calculation in fb after cuts using cut parameter Ej, = 200 
GeV. We list the hard scattering px {Pt^) ranges over which the background processes were evalu- 
ated, and then the backgrounds from various SM processes. The upper bounds quoted correspond 

to the one event level. We take jrif = 170 GcV. 
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FIGURES 



FIG. 1. Contour plots of squark and gluino masses in the uiq - mi/2 plane of the mini- 
mal SUGRA model. Frames arc shown for a) tan/3 = 2, /x < 0, &) tan/3 = 2, fi > 0, c) 
tan/3 = 10, fj, < 0, and d) tan/3 =10, fi > 0. We take nit = 170 GcV and = 0- The bricked 
regions are excluded by theoretical constraints discussed in Paper I, while the shaded regions are 
excluded by experiment. 

FIG. 2. Same as Fig. 1, except we plot contours of lightest chargino mass and contours of right 
slepton mass. Also shown by dotted contours are the kinematic limit for the neutralino spoiler 
decay modes, above which the decays Z2 ZiH^ or Z2 ZiZ are kinematically allowed. 

FIG. 3. SM backgrounds to various SUSY search event topologies in fb, after cuts, but as a 
function of the cut parameter defined in the text. We show frames for a) 1£ + jets events, b) 
OS dilepton + jets events, c) SS dilepton + jets events and d) 3£ + jets events. 

FIG. 4. SUSY signal cross sections for six SUSY cases listed in the text, and total SM back- 
ground in fb, after cuts, as a function of the Ej^ parameter, for the same event topologies as in 
Fig. 3. 

FIG. 5. Contours of cross section (in fb) after cuts described in the text for 11 + jets + I^j, 
events. The solid contours have E^ = 100 GeV, while the dashed contours are for 1,2 and 4 fb 
cross sections with E^ = 400 GeV, from which the maximum reach is derived. The frames are for 
the same SUGRA parameter choices as in Fig. 1. 

FIG. 6. Charge asymmetry Ac defined in the text of the isolated lepton in 1^+jets + ^y events. 
We have fixed ^0 = and tan/3 = 2. 

FIG. 7. The mean jet multiplicity {rij) in 1^ + jets + ^j, events for the same cases as in Fig. 6. 
We have ^0 = and tan/3 = 2. 

FIG. 8. The mean tagged 6-hadron multiplicity {tlb) in l£+jets + -9y events for the same cases 
as in Fig. 6. We have fixed = and tan/3 = 2. The tagging requirements are described in the 
text. 

FIG. 9. Contours of cross section (in fb) after cuts for OS dilepton +jets + ^-j, events. The solid 
contours have E!f = 100 GeV, while the dashed contours are for 1,2 and 4 fb cross sections with 
E^ = 200 GeV, from which the maximum reach is derived. The frames are for the same SUGRA 
parameter choices as in Fig. 1. 
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FIG. 10. Flavor asymmetry (Ap) defined in the text for the OS dilepton + jets + I^rp event 
sample for the same parameters as in Fig. |. We use = 200 GeV (denoted by diamonds) 
except when mi/2 = 100 GeV for which we use = 100 GeV (denoted by squares). The hollow 
(filled) symbols denote 0.2< Ap < 0.5 {Ap > 0.5), while crosses show the points sampled for which 
Ap < 0.2, which is consistent with zero in our simulation. 

FIG. 11. Same as Fig. 9, except for SS dilepton + jets + events. 

FIG. 12. Charge asymmetry Ac, defined in the text, of SS dileptons in i^l'^+jets+^rp events. 
We have fixed Aq = and tan /? = 2. 

FIG. 13. Same as Fig. 9, except for 3£ + jets + I^rp events. 

FIG. 14. An illustrative example showing the variation in cross section after cuts versus the 
SUGRA parameter Aq, for 1£, SS, OS and 3i + jets + I^j^ events. Other SUGRA parameters are 
listed in the figure. We take E!f. = 100 GeV. 

FIG. 15. Regions of the uiq vs. mi/2 plane where clean (central-jet vetoed) isolated tri-lepton 
events are likely to be observable at the LHC, assuming 10fb~^ of integrated luminosity. The 
frames are the same as in Fig. 1, except for the mi/2 scale limits. The filled boxed correspond 
to a lOcj effect above background, open boxes to a 5a effect, and crosses correspond to sampled 
points which were not observable with 10 fb^^. In addition, the kinematic boundary for various Z2 
two-body decays are shown. 

FIG. 16. Regions of the mo vs. mi/2 plane where clean (central-jet vetoed) isolated dilepton 
events (usually from slepton pair production) are likely to be visible using hard slepton cuts de- 
scribed in the text, assuming 10fb~^ of integrated luminosity. The frames are the same as in Fig. 1, 
except for the scale limits. The various symbols correspond to the cross section levels after cuts 
listed on the figure. The estimated SM background level is 0.07 fb. The solid contours correspond 
to m^^ = 70, 100, 150, 200 and 250 GeV, increasing from the lower left. 

FIG. 17. Regions of the mo vs. mi/2 plane where clean (central-jet vetoed) isolated dilepton 
events (usually from slepton pair production) are likely to be visible using soft slepton cuts, as- 
suming lOfb"^ of integrated luminosity. The frames are the same as in Fig. 1, except for the 
scale limits. The various symbols correspond to the cross sections after cuts at the < 3a, (3 - 5)a, 
(5 - 10)<T, (10 - 20)cT and > 20cr levels. The sohd contours correspond to m^-^ = 70, 100, 150, 200 
and 250 GeV, increasing from the lower left. 
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FIG. 18. A summary of the LHC reach (assuming 10fb~ of integrated luminosity) in the 
tuq vs. mi/2 plane for the four cases of Fig. 1 via the various multilepton channels discussed in this 
paper. The dashed-dotted curves show the maximal LHC reach (obtained for some choice of Ej,) 
for the li, SS, OS and 3£+jets + ^/p signals. Also shown is the reach via the complementary clean 
dilepton (marked £) and clean trilepton (marked W1Z2) channels. The boundary of the parameter 
plane that can be probed via va\i\i\]Qi+Ifrp events (with no isolated leptons, denoted by I^j,) as 
obtained in Ref. ||ll| shown for comparison as the dashed curve. 
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